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We revisit the problem of wavefunction statistics at the Anderson metal-insulator transition (MIT) 
of non-interacting electrons in d > 2 spatial dimensions. At the transition, the complex spatial 
structure of the critical wavefunctions is reflected in the non-linear behavior of the multifractal 
spectrum of generalized inverse participation ratios (IPRs). Beyond the crossover from narrow 
to broad IPR statistics, which always occurs for sufficiently large moments of the wavefunction 
amplitude, the spectrum obtained from a typical wavefunction associated with a particular disorder 
realization differs markedly from that obtained from the disorder- averaged IPRs. This phenomenon 
is known as the termination of the multifractal spectrum. We provide a field theoretical derivation 
for the termination of the typical multifractal spectrum, by combining the non-linear sigma model 
framework, conventionally used to access the MIT in d = 2 + e dimensions, with a functional 
renormalization group (FRG) technique. The FRG method deployed here was originally pioneered 
to study the properties of the two-dimensional (2D) random phase XY model [D. Carpentier and 
P. Le Doussal, Nucl. Phys. B 588, 565 (2000)]. The same method was used to demonstrate the 
termination of the multifractal spectrum in the very special problem of 2D Dirac fermions subject 
to a random Abelian vector potential. Our result shows that the typical multifractal wavefunction 
spectrum and its termination can be obtained at a generic Anderson localization transition in d > 2, 
within the standard field theoretical framework of the non-linear sigma model, when combined with 
the FRG. 



I. INTRODUCTION 

Quantum interference induced by multiple elastic im- 
purity scattering can produce very complex spatial fluc- 
tuations in electronic wavefunctions. The statistics of 
these fluctuations may be used to distinguish different 
regimes of qualitative wavefunction behavior, e.g. local- 
ized versus extended. Of particular interest are the wave- 
function statistics at a delocalization transition, such as 
the Anderson metal-insulator transition (MIT) 1 - at the 
mobility edge in three spatial dimensions'' or the in- 
teger quantum Hall plateau (IQHP) transition in two 
dimensions'" Here, the spatial structure of the criti- 
cal wavefunctions is known not to be characterized by 
just a single (or a few) independent exponcnt(s), but 
by an infinite set thereof ('multifractality'). More pre- 
cisely, wavefunction statistics are encoded through the 
r(q) spectrum, or its Legendre transform, the singularity 
spectrum /(a)""" 

The r(g) spectrum is defined via the (generalized) in- 
verse participation ratio (IPR)' given by 

P g ( Si )= [ d d r\^(r)\ 2 ", P g ~L- T M, (1.1) 

JL d 

where d is the spatial dimensionality of the system, L d 
denotes the system volume, and |-0i(r)| 2 is the probabil- 
ity density of a normalized eigenstate wavefunction ipi(r) 
with energy £j, evaluated at the point r. For eigenener- 
gies e lying within a band of extended plane wave states, 



r(q) = d(q — 1), while exponentially localized states yield 
r(q) ~ for L ^> £, with £ the localization length. Multi- 
fractal behavior refers to non-linear g-dependence of the 
r(q) spectrum, and occurs, e.g., at the mobility edge 
e — e c in a disordered three-dimensional (3D) system 
of non-interacting electrons' The singularity spectrum 
f(a) is related to the r(q) spectrum through the Legen- 
dre transformation, 

f(a)=qa-T(q), ^ = a. (1.2) 
aq 

The set of points at which an eigenfunction takes the 
value |V"(r)| 2 ~ L~ a is distributed according to the 
weig ht in this sense, the singularity spec- 

trum characterizes the interwoven fractal measures of the 
sample associated with differently-scaling components 
of wavefunction intensity. The wavefunction statistics 
have been studied experimentally using thin microwave 
cavities p» a very broad distribution of the wavefunction 
intensity, indicative of multifractal behavior, was indeed 
observed. 2 ' 6 ' 12 

The multifractal spectrum [r(q) or f(a)] at a delocal- 
ization critical point is universal, and thus serves as a 
"fingerprint" of the spatial structure of wavefunctions. 
Spectra have been computed numerically at myriad de- 
localization transitions occu ring in variou s spat ial di- 
mensions; see e.g. Refs. lIlilililliili^IilliiHIlil. 

In particular, extensive numerical studies of the IQHP 
transitio n 5 ' 14 ' 16 employing different microscopic models 
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have convincingly established the universality of the en- 
tire f(a) spectrum. Recent work includes that of Refs. 
[lH and which aim in part at decrypting the critical 
(conformal field) theory describing the plateau transition. 

To compute the entire multifractal spectrum analyti- 
cally is, however, a very difficult task in generic systems. 
This is even more so because it is a non-analytic function 
of q or a. As emphasized in Rcfs.[22], [23|, andy, this non- 
analyticity is related to the fact that the r(q) and f(a) 
spectra arc defined for a typical representative wavefunc- 
tion, drawn in principle from a system in a single, fixed 
realization of the static disorder. On the contrary, an- 
alytical methods (i.e., those based upon field theories) 
are best suited for calculating quenched averaged quanti- 
ties. To be precise, we define, following Ref. [23|, two sets 
of multifractal statistics in terms of the 1PR defined in 
Eq. (fTTTl) : 



In this paper we calculate the typical multifractal spec- 
trum at the Anderson MIT in the unitary^ [broken time- 
reversal] symmetry class of disordered, normal metals, 
in d > 2. The spectrum f(q) associated to the aver- 
aged IPR, evaluated at the metal-insulator transition in 
d = 2 + e, was obtained long ago 8 ' 32 ' 33 ' 34 via standard 
perturbative renormalization group (RG). The form of 
the typical r(q) has been argued before only on heuris- 
tic grounds We compute here for the first time the 
typical spectrum directly, using an (analytical) func- 
tional renormalization group (FRG) schem o 25 i 26 ' 27 pre- 
viously employed in the study of wavefunctions statistics 
in a special class of disordered Dirac fermion models in 

2D 22 - 24 - 2 5-26.27.35.36.37.38.39.40.41 



A. Average vs. typical spectra and termination 



r(q) 
HQ) 



dlnP q 
dlnL 
d \sr~P~ q 
d\nL 



(1.3a) 
(1.3b) 



In this equation, the overbar 7TT represents an average 
over realizations of the quenched disorder. The typical 
r(q) spectrum in Eq. (JTTTJ) obtains from the log of the 
IPR for a representative wavefunction; since the latter 
quantity is expected to be self-averaging at the dereal- 
ization transition , 15 i 23 we may introduce an additional, 
though redundant ensemble average over disorder real- 
izations, as in Eq. (|1.3a|) . We have also defined f{q) in 
Eq. (|1.3bl> . which obt ains from the average of the IPR 
itself. The averaged IPR can be encoded through the 
moments of the local density of states (LDOS) operator 
in an effective low-energy field theory (see SecHU below); 
then, the scaling dimensions of the LDOS moment oper- 
ators directly determine f(q). No such effectively local 
construction exists for the typical spectrum r(g), and in 
fact "non-local" (or more precisely, "multilocal" ) correla- 
tions play an essential rol o 24 i 25 ' 26 i 27 in the "termination" 
(defined below) of the typical r(q), as we show in this 
paper. 

For not too large one expects that 



T(q) = f(q), 



(1.4) 



which is the case when the IPR P q represents a self- 
averaging quantity [see Subsection QUI C|) for a re- 
view]. At sufficiently large \q\, however, P q becomes 
broadly distributed ;2 i 12 ' 15 i 22 ' 23 i 28 and the corresponding 
f(q) spectrum, dominated now by "rare events" in- 
duced by the disorder averaging procedure, deviates from 
r(q).— While f(q) is always easier to evaluate analyti- 
cally, it is r(q) that is most easily obtained from a rep- 
resentative wavefunction in numerics.—^ By comparison, 
the average f(q) and f(a) spectra were computed only 
recently via numerics at the IQHP 1 ^ and Anderson 1 ^ 
transitions^ 



In the field theory description of Anderson localiza- 
tion [especially the non-linear sigma model (NLtrM) 
formulation^^ 2 - reviewed in Sec. |TT] the exponent f(q), 
q £ N of the averaged IPR can be read off from the scaling 
dimensions x* and x\ of local composite operators O q (r) 
and 0i (r), which represent the q th and 1 st moments of 
the local density of states (LDOS), respectively i 8 ' 43 i 44 



f(q)=d(q-l)+x*-qxl. 



(1.5) 



(See Sec. [IT] for details.) For example, at the Anderson 
metal-insulator transition in d = 2 + e dimensions in the 
unitary symmetry class, one obtain a 32 ' 33 i 34 



* 
x 1 



-~q{q-l)+0[e 2 q 2 {q~l) 2 } 



e/2 + O e 



: 5/2 



(1.6a) 
(1.6b) 

(1.6c) 



We can define a corresponding average singularity spec- 
trum via 

f{a) = qa-T{q), — — = a, 



dq 

d - h ( a - "o) 2 + O \fe(a-oiQ) 



where 



h 



1 

4S 



0(e) : 



q = d + H + O ( e 5/2 



(1.7) 

(1.8a) 
(1.8b) 



The corrections to [O (■ • • ) terms in] Eqs. (|1.6a[) . (|1.6cl) . 
and (fl~71) - (|1.8bj) obtain at the fourth loop order— (or be- 
yond) in the epsilon expansion. By contrast, Eq. (|1.6b[) is 
exact, and is equivalent to the statement that the average 
(global) density of states is non-critical at the MIT in the 
unitary symmetry class i 45 ' 46 i 47 In the present paper, we 
work only to the lowest non-trivial order in the expan- 
sion parameter y/e. The consistency of the e-cxpansion 
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in dealing with high moments of the LDOS operator is 
demonstrated in Sec. IIV1 Results similar to Eqs. (|1.6aD 
(|1.6c[) were first computed for the time-reversal invari- 
ant orthogonal^^ symmetry class. The so-obtained f(a) 
spectrum is consistent with large-scale numerics^ 

If one were to reconstruct the probability distribution 
of the wavefunction amplitudes from the average spec- 
tra [Eqs. (|1.5[) and (|1.7p ]. a quadratic a-dependence of 
/(a) implies log-normal asymptotics of the distribution 
function . 43 i 44 ' 50 i 51 ' 52 The precursor of this broad distribu- 
tion is already visible at the crossover from the ballistic 
to diffusive regime, where wavefunctions start to show 
(weak) Anderson localization.— In this "pre-localizcd" 
regime, renormalization group studies of an (extended) 
NLcrM ) 43 ! 44 as well as semi-classical analyses of the su- 
persymmetric (SUSY) NLctM^^ 1 -^ predict that the dis- 
tribution of the wavefunction amplitudes starts to devi- 
ate from the Gaussian, developing a log-normal tail£^ 
As it obtains from the f{q) spectrum associated with the 
average of the IPR ) 43 ' 44 this tail reflects the influence of 
rare realizations of the disorder and so-called "anoma- 
lously localized states." 50 i 51 i 52 Even though the tail of 
the distribution is still small, describing rare events in 
the mesoscopic regime, it is responsible for anomalous 
current relaxation, which is slower than expected from 
the Drude formula. 44 ' 50 

For small y/e (i.e., weak disorder) and q not too large, 
one might be inclined to expect that the results for f(q) 
and /(a) in Eqs. (|1.5p ~ (|1.8b[) should not differ substan- 
tially from r(q) and /(a), respectively. However, the 
range of applicability of Eq. (jl.7|) to the typical f(a) is 
limited to a- < a < a+, where, to lowest order 



a± = (Vd ± VZ 



(1.9) 



so that f(a±) = 0. For a > ct + and a < a_, the average 
singularity spectrum /(a) becomes negative, which does 
not make sense if it is interpreted for a typical wavefunc- 
tion [see the discussion following Eq. (|1.2[) . above]. These 
thresholds define the critical values qf of q for the r(q) 
spectrum through 



± _ df(a T ) 




±<?c 



(1.10) 



For q > <7+ , q < q~ , the typical spectrum r(g) deviates 
completely from the average f(g), given by Eqs. ()1.5I) 
(|1.6cl) to lowest order in the epsilon expansion. Indeed, 
it can be rigorously proved 4 - that the r(q) spectrum (as 
defined for a typical wavefunction) must be a monotoni- 
cally increasing function of q; by comparison, the average 
spectrum f(q) in Eqs. (|1.5p - (|1.6c[) is monotonically de- 
creasing for q > (d + S)/2S. For q > q^ , q < q~ , the rare 
maxima (minima) of the wavefunction amplitude domi- 
nate the IPR [Eq. (jl.ip ], as computed for a representa- 
tive wavefunction in a fixed disorder realization. In this 



regime, the associated r(q) is linear in q. By contrast, 
Eq. (|1.4p holds for q~ < q < q+ . In Fig. [TJ we plot the av- 
erage spectrum f (q) as given by Eqs. (|1 .5|) (|1 .6c[) . as well 
as our final result for the typical spectrum r(g), which 
we obtain in Sec. IIIII of this paper [Eq. (|3.24p . below]. 

We say that the multifractal behavior of the typical 
r(q) spectrum "terminates" at q = q^r . This result in 
turn implies that the singularity spectrum f(a) must 
also suffer "termination," i.e., vanish outside of the range 
bounded by a±. The paramount distinction between 
typical vs. average spectra is therefore summarized as 
follows: the termination of r(q) and /(a) reflects the 
dominance of rare amplitude extrema occuring in a rep- 
resentative wavefunction computed for a particular con- 
figuration of the disorder, whereas the deviation of f(q) 
and f(a) from the former reflects the influence of rare 
disorder realizations that enter into the averaged IPR, 
P~ q [Eq. (OB]. 



B. Operator product expansion and the functional 
renormalization group 

The dimension x* (q G {1,2,...}) in Eqs. (|1.5|l and 
(|1.6ap describes the scaling of the disorder-averaged q th 
LDOS moment at criticality, represented by the operator 
O g (r). In order to extract the evolution of the typical 
value of an LDOS moment, we require a scaling equation 
for its entire probability distribution: a functional RG 
approach will turn out to be necessary. We will demon- 
strate that the scaling of the typical LDOS moments de- 
termines the r(q) spectrum. 

A key technical difference distinguishing the calcula- 
tion of the typical r(q) spectrum from that of the average 
f (q) spectrum is that different LDOS moments couple to 
each other along the FRG flow. This coupling among 
the moments is encoded in the operator product expan- 
sion (OPE) of the scaling operators at the derealization 
critical point, 



O q (r)O q ,(r') 



|r - v '\ x « + '<-- x * q+q - q+q ' 



(1.11) 

Whenever the OPE coefficient C q+ q q ^ 0, lower moments 
generate higher ones upon the RG transformation. The 
{x*} satisfy the convexity relation^ 3 - 



-q+q' 



< x„ 



x* n , < 



(1.12) 



for q, q' > 1. Since the {x*} are negative here, Eq. (|1.12p 
indicates that higher moments are much more relevant, 
and hence we are forced to retain all mutually coupled 
moments in the theory, without being able to resort to 
truncation. The FRG allows us to organize and track the 
entire infinite tower of LDOS moment operators. The 
non-zero OPE coefficient C«+? leads to a non-linearity 
within the FRG; the unbounded broadening suggested by 
the g-dependence of x* (reflecting the ever more relevant 
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FIG. 1: Sketch of the multifractal spectra at the unitary class 
Anderson MIT. In the top panel, the heavy solid line corre- 
sponds to the average spectrum f(q), defined by Eq. (|1.3b[) 
in the text, as obtained at the lowest non-trivial order in 
the e-expansion [Eqs. (|1.5j) - (|1.6c[) ] ,22i22i2£ For comparison, the 
heavy dashed line in the same plot is the linear spectrum for a 
plane wave state, f(q) = r(q) = d(q—l). In the bottom panel, 
the solid and dashed heavy line segments represent the typical 
spectrum r(q), defined by Eq. (|1.3a[) . as obtained in this pa- 
per via the functional renormalization group [see Eq. Q3.24)) ]. 
For <k~ < q < q^ (solid segment of the curve in the bottom 
panel), the typical and average spectra coincide. By contrast, 
the typical spectrum is linear for q < q~ , q > q^ (beyond 
"termination"), as depicted by the dashed curve segments in 
the bottom panel. The two curves in the top panel are ren- 
dered as faint gray lines in the bottom, for comparison. The 
inset in the top (bottom) panel depicts the average (typical) 
singularity spectrum at the unitary class MIT corresponding 
to the as-sketched f(q) [r(q)]. 



nature of the corresponding operators, with increasing q) 
is balanced by this non-linearity. For small enough values 
of q, the non-linearity will entirely offset the unbounded 
broadening and render it inconsequential, whereas for 
sufficiently large values of q this will result in the ter- 
mination of the typical r(q) spectrum. 

The mechanism described above is known to be respon- 
sible for the termination of the multifractal spectrum 
in a special (so-called 'chiral')^ symmetry class of 2D 
models, possessing quenched disorder. Carpentier and 
Le Doussa l 25 ' 26 pioneered the use of the FRG technique 
in their study of the random phase XY (gauge glass) 
model. This method was later applied 2 ^ to the problem^ 5 - 
of a 2D massless Dirac fcrmion, subject to a static, ran- 



dom Abelian vector potential. The FRG provided di- 
rect confirmation of the multifractal termination for this 
problem, a result previously conjecture d 22 ' 37 for the vec- 
tor potential model. Later, Mudry et ali^J. extended the 
FRG to a more general 2D disordered Dirac model be- 
longing to the symmetry class BDl (chiral orthogonal 
symmetry class). [We have adopted the nomenclature for 
quantum disorder classes employed in Ref. [48l ] In these 
works, the FRG equation constructed from the set of op- 
erator scaling dimensions x* and OPE coefficients C^ + q f 
[Eq. (jl.lip ] takes the form of the so-called Kolmogorov- 
Petrovsky-Piscounov (KPP) equation^ which describes 
non-linear diffusion in one dimension. It is the non-trivial 
behavior of the long-time asymptotics of the solution to 
the KPP equation that is responsible for the termination. 
We will show that the same equation arises in the gen- 
eral context of the typical r(q) spectrum in the unitary 
symmetry class, at the Anderson MIT critical point in 
d = 2 + e (with obvious extensions to additional symme- 
try classes). 



C. Outline 

Using the framework of the fcrmionic replica (compact) 
NLerM approach ) 1 ' 42 ! 55 we compute the OPE coefficient 
Cq + q ? at the critical point in d = 2 + e for the unitary 
class. Combining this result with the scaling dimensions 
given by Eqs. (|I.6a|) - (|I.6cp . we formulate the functional 
renormalization group for the tower of LDOS moment 
operators that enters into the computation of the typ- 
ical r(q) spectrum. Then we use the FRG to demon- 
strate that the same mechanism active in the 2D Dirac 
models , 25 ' 26 ' 27 discussed above, leads to the termination 
of the multifractal spectrum at the MIT. We obtain the 
r(q) spectrum for a typical wavefunction, which agrees 
with previous heuristic arguments . 15 ' 23 

The rest of the paper is organized as follows: 

In Sec. HIl we review the connection between the IPR 
and the local density of states, and we introduce a gen- 
erating function that will be used to determine the typ- 
ical t(c/) spectrum. We then establish conventions for 
the fermionic replica NLerM, and identify the compos- 
ite operators that represent moments of the local density 
of states in the low-energy field theory. In Sec. IIIII we 
use the operator product expansion (OPE) of the LDOS 
moment operators at the MIT as input into the FRG, 
which then allows us to compute the scaling behavior of 
the generating function introduced in Sec. HH We thereby 
obtain the typical r(q) spectrum. We discuss our results 
and draw conclusions in Sec. IIVI 

The derivation of the OPE of the operators {O q (r)} 
representing the LDOS moments, which constitutes the 
technical field theoretic content of this work, has been rel- 
egated to Sec. El In this Section, we rederive the anoma- 
lous scaling dimensions of the LDOS moment operators, 
and we compute the required OPE coefficient between 
properly normalized versions of these. The results ob- 
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tained are invoked as needed in the earlier Sec. IIII1 so 
the reader less interested in calculational details may skip 
Sec. W\ entirely. 



II. DEFINITIONS AND MODEL 

Extracting multifractality from the LDOS 
typical spectra 



Consider the local density of states (LDOS), defined 



a.s 



v(e,r) = — tenG R (e;r,r) 

= Y,S(e~s l )\Mr)\ 2 , 



(2.1) 



where the retarded Green's function is given by 

^W^(r') 



e — Si + in 



(2.2) 



with T) — > + . 

On the metallic side of the derealization transition, 
we cannot relate P q {s), defined in terms of a single wave- 
function by Eq. (jl.ip , directly to the LDOS4 In order to 
use the field theory approach, we require that the LDOS 
constitute a smooth, well-defined function of energy in 
a closed, finite-size system; this necessitates the reten- 
tion of the finite energy level broadening n > A, where 
A is the global level spacing. (Although a formal de- 
vice in this context, the broadening may be attributed 
to, e.g., inelastic relaxation processes neglected in the 
non-interacting, single particle approach.) 

We define^ 



1 



/ d d Yvi(e,v) 
L<*(9-i)p(<j)(e) = [Jd d Tv(e,T)]"' 



(2.3) 



The quantity p^ (e) denotes the participation ratio, 
which receives contributions from states with energies re- 
siding in a window of width r\ about e. On the metal- 
lic side of the transition, the right-hand side (RHS) of 
Eq. ((23)) should scale identically as Eq. (fTTTj) £ 

It will prove useful to introduce the moment generating 
function for the g th power of the LDOS, 



F q (£;L) = (exp 



J d d TV q {e,r) 



(2.4) 



where the angle brackets (■ • ■ ) denote a suitable ensemble 
average over realizations of the quenched disorder; L is 
the linear system size. Using the identity 



hi 6= -f 



(2.5) 



and replacing P q with the RHS of Eq. (f2~3)) in Eq. (TQaj) . 
the typical multifractal spectrum exponent r{q) may be 
written as 



r{q) 



dhiL 



[F q (C,L)-qF 1 (C,L)}. (2.6) 



Our goal is to compute the scaling behavior of the mo- 
ment generating function F q (£;L), and thereby obtain 
the typical r(q) spectrum via Eq. (|2.6[) . In closing this 
subsection, we note that the evaluation of Eq. (|2.6p using 
the lowest order cumulant expansion for F q (£;L) recov- 
ers the average f(q) spectrum, Eq. (|1.5|k we will discuss 
this point in detail in Sec. IIIII 



B. NLctM formulation 

We examine in this paper the properties of the mul- 
tifractal spectrum at the Anderson MIT in the unitary 
symmetry class. The critical point itself is accessed via 
the standard perturbative e-expansion in d = 2 + e di- 
mensions, with < e <C 1. Our low-energy, effective field 
theory starting point is the compact replica NLcrM^ i 42 i 55 
defined by the functional integral 



V\ 



where 

S[Q] = Y t J ^ rTr (VQ • VQ) - hi d d rTr (k z Q 

(2.7) 

In this equation, the "temperature" t is inversely pro- 
portional to the dimensionless dc conductance of the dis- 
ordered metal, while the "external field" h serves as an 
infrared regulator, coupling to the local density of states 
(LDOS) operator, as defined below. The symbol Q de- 
notes a In x 2n Hermitian matrix field satisfying 



Q 2 (r)=i 2 „, Trg(r)=0. 
The constant matrix 

= diag 



(2.8) 



(2.9) 



sets the (trivial) saddle-point for the action defined by 
Eq. (|2.7p . The identity in the space of 2n x 2n and n x n 

square matrices is denoted by l2n and I„ in Eqs. (|2.8p 
and (|2.9p . respectively. In these equations, n is propor- 
tional to the number of replicas, with n — > at the end of 
the calculation.^^ The target space of the NLcM is the 
compact coset G(2n)/G(n) x G(n), where G = 0,U, Sp 
for the orthogonal, unitary, and symplectic symmetry 
classes, respectively. In the following, we focus upon 
the unitary universality class, G = U. The field the- 
ory in Eqs. (|2.7p and (|2.8p can be derived^ from a mi- 
croscopic Grassmann path integral describing a system 



G 



of non-interacting fcrmions, lacking timc-rcvcrsal invari- 
ance, averaged over configurations of a Gaussian, white 
noise-correlated random potential. 

We employ l a-ir' coordinates^ on the target manifold, 



Q = 



In 



1/2 



w 



1/2 



(2.10) 



For the unitary class, W(r) — > W a p(r) is an uncon- 
strained, complex- valued matrix, with a, j3 € {1, . . . , n}. 

Non-interacting electrons residing in d > 2 spatial di- 
mensions and subject to quenched disorder possess a dif- 
fusive metallic phase, defined as the presence of extended 
wavefunctions at the Fermi energy, provided that the dis- 
order is sufficiently weak. The disorder strength is quan- 
tified by the "bare" conductance at the scale of the mean 
free path, proportional to l/t in the effective field the- 
ory [Eq. P7P]. In direct analogy with the 0(3)/0(2) 
NLcrM description of classical magnetic ordering ) 56 ' 57 i 58 
the "low temperature" (weak disorder) regime < t < t* 
of the model in Eq. (|2.7p exhibits spontaneous continuous 
symmetry breaking, so that the 'cr' fields (I„ — WW^) 1 / 2 
and (I n — W^W) 1 / 2 , which form the diagonal elements 
of the Q-matrix in the parameterization of Eq. (|2.10l) . 
acquire non-zero expectation values throughout the dif- 
fusive metallic phase. By contrast, the off-diagonal i w' 
fields W and represent small spatial fluctuations with 
vanishing mean in this regime. Here, t = t* > locates 
the MIT in d = 2 + e. 

An unusual aspect of the theory of the MIT transcribed 
in Eq. (|2.7p is the fact that this spontaneous symmetry 
breaking occurs also at the derealization transition itself 
(t = t*), and survives even into the insulating ("high 
temperature") phase (t > t*)M^- In the effective NLcrM 
field theory, the trace of the matrix A z Q(r) represents 
the LDOS v(e,r) [Eq. (|2.ip ] for the disordered electron 
system^ 



v(e, r )~Tr A z Q(r) 



Tr 



1/2 



l n -WW\r) + l n -W^W{v) 



1/2 



(2.11) 



This is the same operator that appears in the action 
Eq. (|2.7[) . where it couples to the external field param- 
eter h. While the character of the typical wavefunc- 
tion changes from extended to localized upon travers- 
ing the mobility edge, as encoded by, e.g., the typical 
multifractal exponent r(q) for q > 2 [Eq. (|2.6p ]. the av- 
erage density of states does not exhibit critical behavior 
across the transition^ The LDOS operator on the RHS 
of Eq. (|2.11[) retains a non-zero expectation value so long 
as the average density of states is non- vanishing; conse- 
quently, the (J-matrix cannot be interpreted as an order 
parameter for the MIT. Technically, this result (an excep- 
tion to Goldstone's theorem)^ obtains from the NLcrM 
only after the replica limit n — > is taken. 



For any non-zero, integral number of replicas n € 
{1,2,...}, the model in Eq. (|2.7[) also possesses a (dif- 
ferent) second order transition at t = t* L > 0, separating 
a low temperature "ferromagnetic" phase (t < t*) from 
the high temperature "paramagnet" (t > i*). In con- 
trast to the replica limit n — > appropriate to the de- 
scription of electronic wavefunction (dc)localization, the 
NLcrM with n > 1 is characterized by a restoration of 
the symmetry at the critical point between the 'cr' (diag- 
onal) and V (off-diagonal) components of the Q-matrix, 
within the parameterization given by Eq. (|2.10[) . This is 
the conventional behavior expected for a classical statis- 
tical mechanics model describing spontaneous continuous 
symmetry breaking in the vicinity of the critical point. 

Let us assume that we are interested only in properties 
of the NLcrM given by Eq. (|2.7p at the critical point, t = 
t* n . Because the symmetry is restored at the transition, 
for non-zero n we are permitted to make the following 
U(2n) "rotation" from A z to A x in Eq. (|2~TT|) : 



^(e,r)~Tr A z Q(r 



Tr 



A x Q(r) 



Tr\w(r) + W\r)], (2.12) 



where denotes the block Pauli matrix generalizing 
Eq. (|2.9p . in the standard basis. 

In the technical field theoretic portion of this paper, 
Sec.El we employ the NLcrM defined by Eqs. ([2~7| - (|230l) 
to extract the properties of the LDOS operator and its 
moments. Our strategy is to work, as usual, at fixed, 
integral n > 1 throughout the intermediate stages of our 
computations. At the critical point in d = 2 + e, we 
are then free to employ the LDOS representation given 
by the RHS of Eq. lf2TT2"]) , Only at the end of our work 
will we perform the required analytic continuation n — > 
(which smoothly deforms i* — * t*), so as to obtain 
(perturbative) results appropriate to the MIT. 



C. LDOS moments as composite eigenoperators 

Higher integral moments of the LDOS can be similarly 
represented by local composite operators in the NLcrM. 
The renormalization group (RG) transformation does not 
preserve the form of an operator 



Tr ( W + W ] 



(2.13) 



obtained by taking a power of Eq. (|2.12p . Nevertheless, 
such a structure can be decomposed into invariant eigen- 
operators, each of which possessing an independent scal- 
ing dimension. 

This idea is most easily understood via analogy to 
the simpler 0(3)/0(2) model^^^ to which the field 
theory defined by Eqs. (|2.7p - (|2.10p reduces for the case 
of n = 1 [since U(2)/U(l) x U(l) ~ SU(2)/U(1) ~ 
0(3)/0(2)]. In this NLcrM, the target manifold is sim- 
ply the two-sphere, parameterized by the unconstrained 
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transverse coordinates n± = tt x ± in y , with z-componcnt 
(7 = y/1 — 7r + 7r_. A complete basis of local eigenopera- 
tors with no derivatives is the set of ordinary spherical 
harmonics {Yj jm (7r + , 7r_, a)}. All operators belonging to 
a given irreducible representation of the symmetry group 
possess the same renormalization; therefore, any linear 
combination of spherical harmonics sharing a common I 
value constitutes an eigenoperator. The field coordinates 
tt± arc themselves eigenoperators belonging to I = 1, as 
is the combination 



v = n+ + tt_ oc Y\,—i —Yi. 



(2.14) 



For an arbitrary integer moment of v, one can use angular 
momentum addition to establish the decomposition 



(«- ++ff _) , =E°y ) . 

where the eigenoperators are defined via 



°f= E ^ m Yj, m (n + ,^,a). 



(2.15) 



(2.16) 



m=-j 



with certain coefficients «y„i- For the highest total an- 
gular momentum block j = I in Eq. (|2.15|) . one has 



(2.17) 



since the "highest and lowest weight states" tt 1 + and tt[_ 
are eigenoperators proportional to Yij and Yi t —i, respec- 
tively. 

The coefficients {n^ m } on the RHS of Eq. (gHH) 
are determined entirely by group theory (i.e., are com- 
posed of sums of products of appropriate Clebsch-Gordan 
coefficients)^ up to an overall m- independent normal- 
ization for all operators belonging to a given total an- 
gular momentum block j. This normalization can be 
established via the convention 



Y,.i = A, V 



(2.18) 



In a similar fashion, the operator in Eq. (|2.13[) should 
be decomposed into a sum of terms belonging to dif- 
ferent irreducible representations of the group U(2n). 
Each such term can be further decomposed into a linear 
combination of basis operators with appropriate "mag- 
netic" quantum numbers determined by the transforma- 
tion properties under the subgroup U(n) x U(n). 

It is useful to push this analogy a little further. In 
order to extract the typical r(q) spectrum in the uni- 
tary class model, we need the scaling dimension of the 
most relevant eigenoperator (in the RG sense) contribut- 
ing to each of the p th LDOS moments in Eq. (|2.13[> . 
pg {1,2,.. .}, as well as the operator product expansion 
(OPE) between pairs of such most relevant eigenoper- 
ators. The most relevant eigenoperator contributing to 



the decomposition of Eq. (|2.13|) . for a given fixed p, is 
analogous to the highest (total) angular momentum op- 
erator o\ [S> contributing to the I th moment of (tt + + 7r_ ) 
in Eq. (|2.15P i^ with I = p. [Precise definitions of the 
eigenoperators that we employ in the U(2n)/U(n) x U(n) 
NLctM are given by Eqs. ([2723]) and |2~24|) . below.] In 
the 0(3)/0(2) model, we can effectively trade the op- 
erator of \ which for large I is a complicated sum of 
many terms according to Eq. (|2.16[) . for its lone "lowest 
weight state" component Y t _i [Eqs. (|2~T7)) and (|2~18)) ]. 
Obviously, both operators share the same scaling dimen- 
sion. Moreover, the structure of the OPE between of 

and of, ' follows from that of the product between their 
lowest weight state constituents. Consider the following 
OPE at zero coupling (i = 0): 



1,1' ^W^l+V ' 



(Y 



i,-i 



■){Yi;- 



where we have defined the OPE coefficient 



J+i' = 
-i,i' — 



A 



i+i> 



■■), 
(2.19) 



(2.20) 



The crucial point is that the relative weight of each term 
appearing in the expansion for the eigenoperator of 
[Eq. (|2.16[) ] is entirely fixed by group theory; only the 
overall, ^-dependent normalization is arbitrary. The re- 
quired OPE coefficient in Eq. (|2 . 20[) is then determined 
by just this normalization for the lowest weight state op- 
erators, Eq. (|2.18p . Of course, this argument neglects 
loop corrections, which may modify the value of the OPE 
coefficient given by Eq. (|2.20p . computable systematically 
within the e-expansion. This, however, cannot alter the 
structure of Eq. (|2~ll)]) . 

With the above in mind, we consider the component 



TiW 



(2.21) 



of the LDOS moment in Eq. (|2.13p . As opposed to the 
sphere model discussed above, this pure W power does 
not represent an eigenoperator for n > 1. However, a 
useful subset^ of the RG eigenoperators can be built out 
of p-fold products of '7r' (W a p) field matrix elements: 



(9 



[r) = -W a \ Pl W 



02 



(2.22) 

where (• • • )y means a suitable symmetrization prescribed 
by a Young tableau y. For fixed p, the most relevant 
operator (in the sense of the RG, at the MIT in d = 2 + e) 
is given by the totally antisymmetric Young tableau,^^ 



O 



a ia2 ...a p , s 



sgn(P) [W a ^ p(1) ---W a ^ p(p) ], (2.23) 
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with P a permutation of p symbols; sgn(P) denotes the 
sign of the permutation. Because of the antisymmetriza- 
tion requirement, each distinct operator defined through 
Eq. (|2.23p is identified by any permutation of a complete 
set of indices {cm} satisfying ai ^ 02 ^ . . . ^ a p , and 
similarly for the {Pi}- Indices range from 1 to n, so that 
many different operators can be associated to each inte- 
gral moment of the LDOS, at least for sufficiently large 
n. 

Physically, we would like establish a one-to-one corre- 
spondence between the p th LDOS moment [v(e, r)] p in 
the disordered electron system, and a single, unique op- 
erator O p in the NLctM field theory that represents its 
most relevant component. This can be accomplished by 
tracing over pairs of indices in Eq. (|2.23|) in the following 
fashion: 



0,(r) = J2 



Q-l —1 



a 



n 

E 

=1 



O ' 



V [aia 2 . 



(2.24) 



With this definition, the eigenoperators 



O2 = 



Tr 


[w) 


1 2 

-Tr 


(*>) 




(20 2 



o. 



Tr W 



- 3Tr ( W 2 ) Tr ( W) + 2Tr ( W 3 



etc., are easily recognized as natural deformations of 
the LDOS moments obtained by taking powers of 
Eq. (j2T2Tj) £i. Moreover, we will establish in Sec. [V] that 
the set {Op} closes under the operator product expansion 
(OPE), up to less relevant operators generated on the 
right-hand side of Eq. which we may ignore. This 

is a sufficient condition to apply the functional renormal- 
ization group method. 

In summary, the operators defined by Eqs. (|2.23p or 
(|2.24p constitute the most relevant component (s) of the 
p th moment f tn e LDOS&22£i at the MIT, and hence 
dominate its scaling behavior there. 



D. Augmented NLcrM 

At the metal-insulator critical point, the scaling of 
the average IPR P q [i.e., the multifractal exponent 
f(q), Eq. (|1.3b[) ] can be extracted solely from the scal- 
ing dimensions x* of the local composite operators 



O 



(r) or O p (r), with p e — this is the 

content of Eq. (|1.5p in the Introduction. By contrast, 
the probability distribution functions of the IPR and 
LDOS [reflected by the typical multifractal exponent 
r(q), Eq. (|1.3a[) ] are described by the complicated gener- 
ating function F g (^;L), introduced in Eq. (|2.4p . In the 
low-energy theory, Fi(£;L) can be represented by the 
NLcrM in Eq. (|2.7p with a bare non-zero external field 
parameter ho given by 



Performing a renormalization group transformation upon 
the NLcrM with ho 7^ generically produces higher pow- 
ers of the LDOS operator as new perturbations to the 
action S, so that terms of the form 



5S = -Y m J d d r{Tr [A 2 Q(r)j } 



(2.26) 



for example, will be generated. Here, Y m is a coupling 
constant. The structure in Eq. (|2.26[) is not invariant un- 
der the RG; with further iterations, it will (a) mix with 
other terms sharing the same "engineering" dimension, 
and (b) fuse with other terms and with itself to produce 
new perturbations. Among the flood of structures that 
arise, we will focus only upon the most relevant terms 
that determine the leading scaling behavior for the gen- 
crating function Fi(£; L) of the LDOS and its moments. 

Anticipating the generation of higher moments upon 
renormalization, wc should augment the action in 
Eq. (|2.7| (with ho — 0) by a term of the form 



5S 



00 „ 

Y,Y P / d d rO p (r), (2.27) 
P =i J 



(2.25) 



where the "traced" moment operators O p were defined 
above by Eq. ([2~24]) . 

At tree level, the operators defined by Eqs. (|2.23p and 
(|2.24[) are dimensionless, so that the corresponding cou- 
pling constants {Y p } are strongly relevant perturbations 
to the NLtrM action. As discussed in Sec. II A| they prove 
even more relevant at the non-trivial fixed point (pertur- 
batively accessible Anderson MIT). Moreover, the higher 
moments are more relevant compared to the lower ones 
[Eq. p.I2p ] . The FRG approach tracks the scaling behav- 
ior of this entire tower of operators, and uses this data to 
make non-trivial predictions about observable statistics, 
such as the typical LDOS. Within the FRG framework, 
only two pieces of information are needed: first, the scal- 
ing dimensions of the operators in Eqs. (|2 . 23[) and (12.241) . 

and second, the coefficient C^, 9 for the operator product 
O q <g> O q , O q+ql , as defined by the OPE in Eq. (TTTTj) , 
All quantities are to be evaluated at the MIT in d = 2 + e. 

We use a two-stage approach to the renormalization 
of the "extended" NLcrM [the action Eq. (f2~7|) supple- 
mented with (|2.27[) ]. The idea is to first locate the non- 
trivial metal-insulator fixed point in d = 2 + e, obtained 
via the standard e-expansion by renormalizing the theory 
in Eq. (|2.7p with h — ► 0. (We will use dimensional reg- 
ularization.) We then compute the OPE [Eq. at 
the MIT to the lowest non-trivial order in s/i. Finally, we 
run a 'one-loop' RG calculation at this non-trivial fixed 
point, for the full model defined by Eqs. (flTT)) and (f!T2T|) . 
The required one-loop functional renormalization group 
equation is obtained from the OPE4^ Note that since we 
are interested in LDOS and IPR statistics at the Ander- 
son metal- insulator transition (t — £*), rather than in the 
diffusive metallic phase (t < V), we are required to run 
the FRG at this non-trivial fixed points 
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In order to streamline the presentation, the above- 
described field theory calculations are relegated to the 
last Sec. |V] of this paper. The obtained results required 
for the functional RG are simply invoked as needed in 
the next Sec. IIII1 so that the reader less interested in 
calculational details may avoid Sec. IVl entirely. 



III. FUNCTIONAL RG FOR THE TYPICAL r(q) 
SPECTRUM 



A. From coupled RG to KPP equations 

The typical t(c/) spectrum, defined in Sec. U by 
Eq. p.3a[) . can be extracted from the generating func- 
tion F q (£;L), introduced in Eq. (|2.4p . The relationship 
is expressed by Eq. (|2.6p . In terms of the NLcrM formula- 
tion reviewed in Sees. Ill Bl III Dl F q (£; L) may be encoded 
as 



^ Y pq I d d ~ 



[Op,](r) 



(3.1) 



where q — 1,2,3, and \O pq \ (r) is a "renormalizcd 
and normalized" LDOS moment eigenoperator, defined 
by Eq. (|Q5|) in the technical Sec. [V] of this paper. [Note 
that here pq denotes the product of the integers p and q. 
[O m J (r) is just a normalized version of the LDOS mo- 
ment operator C m (r), defined previously via Eq. (|2.24l) . 
The careful normalization of operators is an important 
technical step required for the accurate computation of 
correlation functions at the MIT, as detailed in Sec. [V] 
In this Section, we merely assert that the proper pro- 
cedure has been implemented.] Eq. (|3.ip generalizes 
Eq. (|2.27p for the case of q > 1: in order to compute 
F q (£; L), one must augment the bare sigma model action 
with the operator tower { [O g ] , [02<j] , [03<?] , ■ ■ since 
through the OPE operators representing lower integral 
LDOS moments generate new ones representing higher 
integral multiples of these. For q > 1, the operators in 
Eq. (|3~T]) form a subset of those in Eq. (j2~2"T)) . The ex- 
pectation (■ • • ) in Eq. (|3.1[) is taken with respect to the 
NLcrM action at the MIT in d = 2 + e, Eq. i[2~7j) . with 
h = and t = t* . The coupling constants Y pq take the 
bare values 



Y pq (l = 0) = -£6 P> 



(3.2) 



Here, I = InL/Lo is the log of the spatial length scale 
L (e.g., the system size), with Lq an arbitrary reference 
scale. 

The simplest approximation to F q (£; L) obtains from 
the lowest order cumulant expansion of Eq. (|3 . 1 [) , evalu- 
ated at L = Lq [i.e. using the bare coupling constants in 
Eq. dO l]: 



Fq(£,;L ) ~ exp 



exp 



(-e^- : ) ■ 



In this equation, x* denotes the negative scaling dimen- 
sion of the operator [Oj (r) at the MIT; for the uni- 
tary class studied here, the result to lowest order in yje 
was given by Eqs. (|1.6a|) and (|1.6c|) , above. Combining 
Eqs. (|3.3p and (|2.6p . we immediately recover Eq. (|1.4p : 
the lowest order cumulant approximation to -Fg(£;Lo) 
equates the typical r(q) spectrum with f(q) [Eqs. (|1.5D 
(| 1 .6c[) ] . associated to the average of the IPR. For suffi- 
ciently large moments with q > q c , where q c was defined 
by Eq. (jl.lOp . this identification invariably breaks down 
(see the discussion in Sec. II Ap ; an accurate computation 
of Eq. (|3.3p then requires the retention of higher order cu- 
mulants. One immediately sees the need for the operator 
product expansion (OPE), as defined by Eq. p. lip : the 
second and higher cumulants involve products of LDOS 
moment operators, integrated over the sample volume. 
When two (or more) such operators approach the same 
spatial position, fusion can occur, in which new, higher 
moment operators are generated through short-distance 
regularization,^^^^^ At the MIT, the negative scal- 
ing dimensions {x*} of the LDOS moment operators 
[Eqs. (|1.6ap - (|1.6cp ] satisfy the convexity relation given by 
Eq. (|1.12p . Therefore, operators corresponding to succes- 
sively higher moments carry ever more negative scaling 
dimensions, contributing ever more strongly to the cu- 
mulant expansion. 

Rather than compute the generating function F q (£; L) 
directly, we will use scaling arguments to extract its 
asymptotic behavior in the large system size limit, 
L/Lq — > oo. In Eq. (|3.ip . the LDOS moment opera- 
tors \O pq \ perturb the action of the critical field theory. 
It is well known that the lowest order RG equations for 
the set of conjugate coupling constants Y pq follow directly 
from the operator product expansion^ 

In Sec. El we demonstrate that (the properly normal- 
ized versions of) the operators defined by Eq. (|2.24p obey 
the OPE given by Eq. (TlTTTj) at the MIT in d = 2 + e. We 
find that the OPE coefficient is given by the "tree level" 
(zero coupling) amplitude 



r<i+i 

U q,q' 



(g + g') 

q\q'\ 



0(e). 



(3.4) 



These results are obtained as Eqs. (|5.33p and (|5.34p in 
Sec. |Vj where we demonstrate that the lowest order 
t* (x y/e (one-loop) correction to the OPE coefficient in 
Eq. (pT4")l vanishes. Using Eqs. (fTTTTj) and (f3"l| . one finds 
the infinite set of RG equations^ 



- (d-x* )Y +— V 

m—l 



(3.3) 



(3.5) 

where Sd is the surface area of the sphere in d dimen- 
sions. Through the OPE, lower moment coupling con- 
stants always generate higher ones; the convexity prop- 
erty in Eq. (|1.12p implies that, for p > p' , a non-zero Y pq 
represents a much more relevant perturbation than Y p > q 
to the critical NLcrM action. Clearly we must retain the 
entire infinite set in our analysis. 
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At first glance, the generation of infinitely many rel- 
evant couplings would seem to imply non-universality: 
there are infinitely many classes of solutions to the RG 
equations (|3.5[) . and hence there are infinitely many ways 
to depart from the RG fixed point representing the MIT. 
This is consistent with the fact that a random critical 
point should be characterized by the entire distribution 
functions of physical quantities, which can become very 
broad. At a derealization critical point, however, the 
multifractal r(q) and /(a) spectra, associated to a typical 
wavefunction in a fixed disorder realization, are both self- 
averagin g 15 ' 23 and universa h 5 ' 14 i 16 ' 18 i 19 We will demon- 
strate that the functional renormalization group (FRG) 
method gives a universal prediction for r(q) and /(a), 
below and above termination (as defined in Sec. II Aj) . at 
the unitary class Anderson MIT in d = 2 + e consistent 
with this picture. 

We can trade the coupled set of ordinary differen- 
tial equations in Eq. (|3.5[) for a single partial differen- 
tial equation (PDE) by defining the auxiliary generating 
functional 

c 00 I —z\p 

G q (z, I) = G q (~z, 1) ee 1 + ^ £ t-J- Ypq (l), (3.6) 

P =i 

where we have introduced the "position coordinate" z. 
G q (z, l) is a "Galilean boost" of G q (z, I), with 

Z = Z + Eql, 1 = 1. (3.7) 

At the unitary class MIT in d = 2 + e, the parameter S 
has the value given by Eq. p.6cp . to one- loop order. 

Using the RG equations (|3.5|) for the coupling 
constants {l^}, and the explicit form of x* q from 

Eqs. (jl.6ap - (|1.6c[) . one can easily show that G q {z,l) 
satisfies the following Kolmogorov-Petrovsky-Piscounov 
(KPP) equation: 

^djG q = D q dlG g + G q (G q -l), (3.8) 

where we have introduced the effective diffusion constant 

D q ee (3.9) 

The same Eq. (|3.8[) was obtained in previous FRG studies 
of 2D disordered systems i 25 i 26 ' 27 



B. Solution to the KPP equation and results 

The KPP equation (|3.8[) describes non-linear diffusion 
phenomena. The positive D q in Eq. (|3.9p reflects the 
diffusion of the distribution function for the inverse par- 
ticipation ratio (IPR), defined by Eq. (jl.ip . For high 
moments, <; > 1. this diffusion constant is very large, in- 
dicating that the IPR becomes broadly distributed in the 
large system size limit; in this regime, the f(q) spectrum 



associated with the average IPR [Eq. (Ob}] is domi- 
nated by rare realizations of the disorder, and loses its 
meaning with respect to the typical wavefunction. The 
non-linear term in Eq. (|3.8p appears because the gen- 
crating functions G q (z,l) and F q (£; L) encode informa- 
tion about the typical r(g) spectrum; this non-linearity 
arises through the OPE between LDOS moment opera- 
tors [Eqs. (jl.lip and (|3.4p ]. As explained in the para- 
graph following Eq. (|3.3p . the OPE is the essential in- 
gredient required in the computation of r(q), which was 
missed in previous treatments 8 - ' 43 ^ 4 of the Anderson MIT 
based upon the NLerM approach. 

Non-linear PDEs are often not analytically solvable, 
but a number of key results are known for the KPP equa- 
tion. We summarize here only those features essential 
to the c omput ation of r(q); for further details, consult 
Refs. [H2i27| and the references therein. For a large 
class of initial conditions which satisfy 

lim G q (z,0) = 1, (3.10a) 
lim G q (z,Q) = Q, (3.10b) 

z — > — 00 

G q {z,l) converges to a stable traveling wave solution 
propagating in the positive z direction, 

G q (z, I — > 00) — > h(z — c q l), (3-11) 

where the constant c q denotes the wavefront velocity. 
The functional form of the traveling wave in Eq. (|3.1ip 
is sensitive to the details of the initial condition at I = 0. 
On the contrary, for an initial wavefront satisfying the 
asymptotic property 

G q (z -> +oo,0) - 1 - Xe~\ (3.12) 

with A a pure number, the velocity c q is universal, de- 
pending only upon the diffusion constant D q , defined in 
the context of the MIT by Eq. (J3J|, above. Note that 
Eq. (|3.12p constrains G q only in the region penetrated 
by the wavefront [Eq. (|3.1ip ] in the limit of large "renor- 
malization time," I — > 00. Remarkably, the wavefront 
velocity is also insensitive to the precise form of the non- 
linear term J-(G q ) = G q (G q — 1) in the KPP equation 
(|3.8p . In fact, the same velocity obtains from KPP for 
any nonlinear forcing function satisfying the constraints 

T(0)=T(l) = 0, J r (G)<0, 

^SU-L ^B = -L (3.13) 
dG ~ dG 

for < G < 1. In this sense, the KPP equation achieves 
a strong version of universality. 

Let us now return to the problem at hand, computing 
the typical r(q) spectrum obtained at the MIT in the 
unitary class for d = 2 + e. The initial condition for the 
KPP Eq. (JSU) implied by Eq. (JS^l) is 

G q (z,0) = l- (3.14) 
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consistent with only Y q non- vanishing. Eq. (|3. 14[) satisfies 
the condition in Eq. (|3.10a[) . having the same form as that 
expressed in Eq. (|3.12|) . In order to satisfy Eq. (|3.10b|) . 
we must bound the amplitude < G q (z — > — oo,0) < 1; 
to that end, wc deform Eq. (|3.2[) as follows: 



Y q (0) 
Y Pq (0) 



Sd 
2d 



p-i 



which leads to 



G q (z,0) ~ exp 



^2d 



(3.15) 



(3.16) 



Crucially, since Eq. (j3~T6|) satisfies Eq. (|3TT2|) . the 
asymptotic traveling wave velocity c q [Eq. (|3.1ip ] de- 
pends only upon the diffusion constant D q , Eq. (|3.9p . 
For t he K PP equation J53J) satisfying (|3.10a|l . (|3.10b[> . 
and (|3.12|) . one finds qualitatively different behavior for 
D q less than or greater than one i 25 ' 26 ' 27 



d(l + D q ), D q <l 



D q > 1 



(3.17) 



Reversing the Galilean boost in Eq. (|3.1ip via Eq. (|3.7p . 
we see that 



G q (z, I — > OO) ^ h(z — C q l) , 

c q = c q - Eq. 



(3.18) 



Let us try to understand the physics implied by 
Eq. (|3.18p . The generating function G q (z,l) was denned 
via Eq. (|3.6p in terms of the infinite tower of coupling 
constants the latter were introduced in the NLctM 

definition of F g (£; L), Eq. (|3.ip . Under a change of length 
scale (e.g. an incremental increase in the sample size), 
each Y pq evolves according to the RG Eq. (|3.5p . If wc 
neglect the non-linear terms in this equation due to the 
OPE, then each non-zero Y pq grows under renormaliza- 

tion according to its own scaling exponent c ; 



(0) 



a full characterization of the system requires the spec- 
ification of the entire set {dpq}, p £ N. As argued in 
Sec. Uand below Eq. (|3.5p . we nevertheless expect that 
a single, well-defined exponent r(q) can be defined for 
the IPR associated to a typical wavefunction, even in 
the limit of relatively "large" q. Eq. (|3.18p implies that, 
through the OPE and the subsequent non-linearity of the 
KPP equation, the functional RG proves this assertion: 
the asymptotic scaling of G q (z, I) involves a single num- 
ber, the velocity c q given by Eqs. (|3.17p and (|3.18p . which 
we can think of as a "typical" scaling exponent. 

Since G q {z, I) tracks the scaling of coupling constants, 
we infer that the set {^g} "fuses" into a single, typical 
coupling Y q yp , up to less relevant perturbations to the 
critical NLerM fixed point; we can then define an associ- 



ated typical anomalous dimension 



x q p — d c q 



-Zq(q-l), 



1 < q <q c 



q (Vd - sgn(q)y/z) , q > q, 



(3.19) 



where we have used Eq. (|3.9p . In this equation, the crit- 
ical value q c was defined in the Introduction, Eq. (jl.lOp . 

We infer from Eq. (|3.19p that F q (£;L) acquires the 
following asymptotic form: 



F q (£; L — > oo) ~ < exp 



Vtyp 

9 



d d r 



exp 



Y q tvp L 



(3.20) 



As in Eq. (|3.3p . we have evaluated F q (^; L) in the lowest 
order cumulant expansion; the crucial difference between 
Eqs. (|3.3p and Eq. (|3.20p resides in the implied order 
of operations. To obtain the final result in Eq. (|3.20p . 
we first coarse grain the system, say by intcgrating-out 
short wavelength degrees of freedom (in a Wilsonian pic- 
ture). The coarse graining generates higher order cou- 
plings {Y pq }, p > 1, through the non-linear RG Eq. (|3.5p . 
In the large system size limit I = \n(L/ Lq) — > oo, a single, 
well-defined typical coupling F g typ emerges, associated to 
a new local operator [Og YP ] (r) , whose scaling dimension 
is given by Eq. (|3.19p . Finally, we evaluate F q (£;L — ► oo) 
to lowest order in the cumulant expansion, which gives 
Eq. (|3.20p . This is expected to be a correct represen- 
tation of the the asymptotic scaling limit, because the 
functional RG has already built all of the most rele- 
vant operator "fusions" into the definition of [C^ yp ](r). 
The emergence of the associated Y^ yp and a;* yp has been 
proven above using the properties of the KPP equation, 

Eq. dm 

Finally, we extract the typical r(g) spectrum. As ob- 
tained in the limit of large, but finite renormalization, the 
typical coupling Y q yp should have an analytic expansion 
in powers of the parameter £ [c.f. Eqs. (|3.2p and (|3.15j) ]. 
Up to an irrelevant rcscaling, we may write 



n typ = -e - E W 



(3.21) 



m=2 



Combining Eqs. (|2.6p and (|3.20p . we obtain 
_d r°° d£ 



r(q) 



qe 



Y, W L 



(3.22) 

In the limit L —> oo, we may neglect all but the first 



term in Eq. (|3.2ip . since d 



r typ 



> for all q > 1, 



provided S < 4c?. This condition is always satisfied in 
the perturbatively ac cess ible regime, < t « 1, where 
the parameter 3 = \/e/2+ O (e) < 1 [Eq. fOc|) ]. Then 
we obtain, using Eq. (|2.5p 



r(q) =d(q - 1) + 



r typ 



typ 

qx{ ■ 



(3.23) 
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Eq. (|3.23|) for the typical r(q) should be compared to 
Eq. p.5p for f(q). In the perturbative regime e <C 1, 
we h ave a :'/ 13 = [Eq. (OS])]. Combining Eqs. ([3l9]l 
and (|3.23p . we arrive at our final result, the typical r(q) 
spectrum given by 



d(q-l)fl--^J, \q\<q c 
T(q) ~ { d(l-^Yq, \q\>q c 



(3.24) 



where q c = Wd/B. [Eq. (|1.10p ] . In this equation, we have 
extended q from the positive integers to the entire real 
line. The average f(q) and typical r(g) spectra are re- 
spectively sketched in the top and bottom panels of Fig.[T] 
in Sec. U 

The singularity spectrum f(a) was introduced in 
Eq. (|1.2|) . The f(a) corresponding to the typical r(q) 
in Eq. (|3"^4]) is 



/(«) 



(a-d- E) 2 
45 

q 2 (a + — a)(a — a_) 
Ad : 



lo, 



a_ < a < a + 

a < a_, a + < a 
(3.25) 



The spectral cutoffs a± were defined by Eq. (|1.9p . As 
expected, f{a) associated to the typical wavefunction is 
never negative, as discussed in Sec. HI see also the top and 
bottom panel insets in Fig.Q] Eqs. (|3~2"4"|) and (j3~2"5)) hold 
to the lowest non-trivial order in ^fe. The consistency of 
the restriction to only the lowest order contributions in 
the e expansion is demonstrated in Sec. IIVI 

An alternative representation of multifractality invokes 
the "generalized dimension" D q , defined via 



r(q) = (q- l)D q 



(3.26) 



Spectral termination [Eqs. (|3.24|) and (|3.25[) ] implies that 
a+ < D q < a_, with the boundary values associated to 
the limits 



lim D„ = azr. 

q— >±oo 



(3.27) 



Numerical computations of D q for the typical wavefunc- 
tion confirm Eq. (|3.27[) ; see, e.g., Refs. [H and [3. 



C. Comparison with other arguments 



Our final results (|3.24[) and (|3.25p ag ree with prev ious 
heuristic arguments given in Refs. l37l [TBI , and l23l. As 
explained below Eq. (|1.2|) . f(a) describes the measure 
£,f( a ) f the set of those points r where the eigenfunction 
ip takes the value \ip(r)\ 2 oc L~ a ^. Hence, the IPR P q 
[Eq. (jl.ip ] can be estimated as an integral 



P,~ / daL- qa+f( - a l 

/(a)>0 



(3.28) 



The integrand takes a maximum value at a saddle point 
value a, which defines the r(q). For \q\ < q c , this saddle 
point is in the integration domain, whereas for |g| > q c 
it is outside of it. In the latter case, the integral is 
dominated by the boundary value of a = a T , where 
/(a T ) = 0. 



IV. DISCUSSION 

We have provided a field theoretical description of the 
termination of the multifractal spectrum r(q), as defined 
for the typical wavefunctions, at the Anderson MIT in 
d = 2 + e for the unitary disordered metal class. The 
essential ingredients of the calculation are evident in the 
formulation of Eqs. (|2.6[) and f)3 ,1ft : these are the infi- 
nite set of properly normalized LDOS moment operators 
{ [O pq j (r)}, p e N, characterized by the negative scaling 
dimensions {x* q } in Eqs. (|1.6ap - (|1.6c[) . Each successive 
higher moment operator with p = {1, 2, . . .} constitutes 
a more strongly relevant perturbation to the critical RG 
fixed point that describes the MIT. Through the OPE 
[Eq. p. lip ] , lower moments always generate higher ones, 
and a consistent treatment of the problem requires that 
the entire infinite hierarchy of LDOS moment operators 
is retained. Ordinarily, the advent of an infinity of rel- 
evant scaling directions should cast serious doubt upon 
the adequacy of single (or few) parameter scaling, at least 
with respect to the investigated critical point; remark- 
ably, the FRG "absorbs" the entire LDOS moment tower, 
and through the (universal) properties of the long-time 
asymptotics of the KPP equation, renders in the end a 
single, universal prediction for the typical r(g). 

Physically, the relevant LDOS moments reflect the fact 
that a random critical point should be characterized by 
the distribution functions of physical quantities, rather 
than their mean, variance, or first few moments. The 
distribution of an observable in the presence of quenched 
disorder can become very broad, due to the influence of 
rare eventSi 43 i 44 ' 50 i 51 ' 52 In principle, wc need the func- 
tional rcnormalization group (FRG) to obtain scaling for 
the entire probability distribution^ 4 - For large q, the IPR, 
defined by Eq. (jl.ip [or of its field theoretic generaliza- 
tion, Eq. (|2.3p ], constitutes such a broadly-distributed 
observable^ i 12 i 15 i 22 i 23 ' 28 By comparison, a universal r(q) 
spectrum for the typical wavefunction obtains because 
the log of the IPR is self-averaging for all g. 15 i 23 

Technically, the FRG method implemented in 
Sec. IIIII is completely analogous to that employed 
previousl y 25 ' 26 ^ 7 in the study of certain special 2D dis- 
ordered field theories, possessing an additional, 'chiral' 
symmetry4£ [See the end of Sec. II B I for a description of 
these chiral models.] As in this prior work, the FRG 
translates the infinite set of coupled flow equations (|3.5p 
into the KPP Eq. (|3.8p for a certain (auxiliary) gener- 
ating function. Through the asymptotic solution of the 
KPP equation in the form of a propagating wavefront, 
the tower of relevant LDOS moment operators combines 
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via multiple OPEs into a single, typical operator (up to 
less relevant perturbations), characterized by the typi- 
cal scaling dimension in Eq. ()3.19j) . The final results for 
r(q) and f(a) [Eqs. (|3.24j) and (|3.25j) ] are obtained via 
the FRG for the unitary universality class, using only 
two inputs, evaluated at the MIT: (i) the scaling dimen- 
sions {x*} (associated to the average operator scaling, 
already known from previous work) 8 ' 32 ' 33 ' 34 and (ii) the 
OPE coefficient C q + q ? , Eq. tf3l} (computed in Sec. ED 
to lowest non-trivial order in ^/e). The former is spe- 
cific to the unitary class, but we have shown that the 
latter takes exactly the same form in the chiral model 
calculations i 25 i 26 i 27 

In treating the unitary class, we have chosen to 
work only to the lowest order in t* oc y/e, i.e. to one 
loop. To this order, the resulting singularity spec- 
trum f(a) given by Eq. (|3.25[) is purely quadratic over 
the region a + < a < a_ (the so-called "parabolic 
approximation") . 4 ' 5 ' 17 ! 18 ' 19 We now discuss the consis- 
tency of working with the functional renormalization 
group to this order in the e-expansion. Corrections to 
the LDOS moment scaling dimensions {x*} are already 
known to four loops ; 33 ' 34 

K = -yj\ g(g - i) - ^V(g - i) 2 + o (V/ 2 ) , 

(4.1) 

where Q{z) denotes the Ricmann zeta function. While the 
FRG method formally retains LDOS moment operators 
{[Cq] (r)} to arbitrarily high orders in q, it is crucial to 
note that the termination of the typical r{q) spectrum 
[Eq. (|3~24"|) ] occurs at the finite value q = q c [Eq. (fTTTTJj) ]; 
to lowest order, 

q 2 c =2^2/~e+0(l). (4.2) 

Evaluating the 4-loop scaling dimension in Eq. (|4.ip at 
q c , we obtain 

x* q=qc = -2 + (2ef 4 +0(e 1 / 2 ). (4.3) 

The one-loop approximation consists of retaining only 
the first term on the right-hand side (RHS) of Eq. (|4.ip . 
as well as the terms written explicitly on the RHS of 
each of Eqs. (|4.2p and (|4.3p . At termination (q = q c ), the 
higher order loop corrections give rise to additional terms 
in Eq. (|4.3p that are down by higher powers of e 1 / 4 , and 
these can be consistently neglected for eCl. 

We cannot resist contemplating, at a very speculative 
level, a naive extrapolation of our one-loop results to 
moderate or even large e. First, note that Eq. (|3.19p 
implies the existence of a non-zero, typical scaling di- 
mension x± p for the first moment of the LDOS, when 
q c < 1: 

x^ = d(l-^j , q c <l. (4.4) 



Using the lowest order result in Eq. (|4.2p , we then define 

(f = e(9c = 1) 

~ 8. (4.5) 

At face value, a x^ p > would imply that the typical 
LDOS vanishes at the MIT. Eq. (|4.5p suggests that this 
becomes possible in the limit of large spatial dimension- 
ality, e > £f- Such a scenario does not contradict rigorous 
result s 45 ' 46 which prove that the average, global DOS re- 
mains uncritical (constant) across the transition for any 
spatial dimension d. Indeed, Bethe lattice computations^ 
exhibit a typical LDOS that vanishes exponentially across 
the transition^ 6 . As the Bethe lattice can be equated with 
the limit of infinite spatial dimensionality, 3 this picture 
in fact appears consistent with Eqs. (|4.4p and (|4.5p in the 
limit e — > oo, where upon x^ p — > oo. 

It has been asserted 3 - that the upper critical dimension 
for Anderson localization occurs only at d — oo, i.e., the 
Bethe lattice case. The naive extrapolation of the results 
of this paper to large e suggests an alternate possibility. 
To motivate the basic underlying idea, we note that for 
e > ep, the typical multifractal spectrum r(q) defined by 
Eqs. (|3~T9|) and (|3~23"]) would take the form 

[ - d (l- ±) , \q\ < q c 
r(q) = ! ^ V I* J (4.6) 

\—{q~\q\)- \q\>Qc 

y q c 

This should be contrasted with Eq. (|3.24p . which as- 
sumed q c > 1 [always the case in the perturbatively ac- 
cessible regime, < e < 1— sec Eq. P~2"]) ]. Eq. (j4"l))) 
shows that r(q) = for all q > q c when q c < 1 
(e > ep). Thus the r(q) spectrum "freezes" for dimen- 
sionalities above the threshold dp = 2 + ep. An analo- 
gous freezing transition has been predicte d 22 ' 25 ' 26 ' 27 ' 37 ' 38 
for the 2D chiral Dirac models discussed at the end 
of Sec. II At in these models, the transition occurs for 
quenched disorder fluctuation strengths larger than some 
threshold value. Unlike the unitary metal class discussed 
here, the chiral model freezing transition has been rig- 
orously derived through strong randomness 3 ^ and FRG 
arguments , 25 ' 26 ' 27 ' 36 which do not rely upon expansion 
in a small parameter. By comparison, Eqs. (|4.4p and 
(|4.6p lie well beyond the perturbatively accessible regime. 
Regardless, the freezing scenario suggests the intriguing 
possibility of a finite dp < oo for Anderson localization 
in the normal metal classes, which could perhaps serve 
as a finite upper critical dimension. 

V. OPERATOR PRODUCT EXPANSION AT 
THE ANDERSON FIXED POINT: 
PERTURB ATIVE CALCULATION 

In this final (technical) section, we provide a deriva- 
tion of the LDOS moment operator algebra required for 
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(a) 



(b) 




FIG. 2: Propagator (a) and lowest-order vertex (b) necessary 
for the one-loop RG, obtained by expanding Eq. (|2. T|) in terms 
of the unconstrained W field, using Eq. (|2.10|l . Associated 
amplitudes are given by Eqs. (|5.1[) and (|5.2|l in the text. 




FIG. 3: Basic renormalization process of a composite opera- 
tor. 



Adopting standard dimensional rcgularization 
conventions ) 57 i 58 i 67 



the functional RG construction in Sec. Mil Using the 
NLcM framework reviewed in Sec. [TTJ we first redcrivc 
the anomalous scaling dimensions of the LDOS moment 
operators introduced in Sec. Ill CI We then turn to the 
perturbative evaluation of the operator product expan- 
sion, as defined by Eq. between properly normal- 
ized versions of these eigenoperators. 



A. Renormalization of the model 

To begin, we consider the renormalization of the 
bare NLctM defined by Eq. ([2~7]) . This calculation is 
standard ^ 42 ' 67 we provide only our conventions neces- 
sary to set up the computation, and the corresponding 
results. Using the parameterization in Eq. (|2.10p . one ob- 
tains the W — > W a p field propagator and vertex shown 
in Fig. [21 The propagator is pictured in Fig. Ufa) as pair 
of counter-directed thin lines, representing physically an 
ambulating electron- hole pair ('diffuson'), and mathe- 
matically the linking of direct and conjugate indices in 
two inequivalent representations of U(n) [since the maxi- 
mum compact subgroup of U(2n) is U(n)xU(n)]. Equiv- 
alently, in the unitary class with this parameterization, 
the field W a p is Wick-contracted only with its adjoint 
W^a, this fact is indicated by the thick arrows in Fig. [2J 
which also encode the direction of momentum flow. The 
amplitude corresponding to the propagator in Fig. Ufa) 
is 



r ft (k)fft«^ 2 (k))^;'^. 



hnt 



n'n 



(5.1) 



where to and ho are bare parameters. In this paper, we 
will only need the lowest order non-linear vertex = QJ4 
[obtained via an expansion of Eq. (|2.7[) in powers of W] ; 
this vertex is pictured in Fig. Hfb) , with the correspond- 
ing amplitude 



2! 

^4 = 7— 

8t 



2(kx -k 3 + k 2 • 
- k 2 • k 3 k 3 



k 4 ) - ki • k 2 
• k 4 k 4 • ki — 2/ioio 



(5.2) 



t = i/i~ £ F t , 

_1 

ho = Z w 2 h, 



(5.3) 



with t and h rcnormalizcd parameters, [i an arbitrary 
inverse-length scale, and Zw the field renormalization of 
the elementary operator W, the RG flow equations are 
given by 



dt 
dlnh 



-et 



-, , din F f ' 
1 ' dint 



1 d In Z\y d In t 
~dT~ ~ + 2 dint dT' 



(5.4) 



e is the spatial dimen- 
— In fi is the logarithm 



In Eqs. (EHJ) and $5A$), d = 2 
sionality of the system, and I <■ 
of the spatial length scale. 

For the compact unitary model with target space 
U(2n)/U(n) x U(n), the one- loop flow equations are 



dt 

Jl 
d In h 



^7 = -et 



47T V ' 



nt 
47 



o (e 



(5.5) 
(5.6) 



These equations possess a critical fixed point at the "tem- 
perature" t* = Ane/n, with h = 0. The Anderson model 
corresponds to the limit n — ► in Eqs. (|5 . 5|) and (|5.6[) ; in 
this case, a non-trivial critical point occurs at two loop 
order ii In our conventions, the two-loop result i a 42 i 67 



dt _ t A 



O t 4 



(5.7) 



valid in the limit n — > 0. The critical value of the inverse 
conductance at the metal- insulator transition in d = 2 + e 
is proportional to t* = Any/2e + O (e). 



B. Composite operator scaling dimensions 

Next, wc turn to the renormalization of the composite 
operators defined by Eq. (|2.23|) . The renormalization of 
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O 



a 1 a 2 ...a p 
p [/3i/3 2 .../3 P ] 



can be determined by considering "matrix 



elements" of that operator with arbitrary configurations 
of p distant, mutually separated adjoint fields {W^ A 7 }. 
[For more general operators built out of products of both 

V (W,Vft) and V (\/i n - WW* ,y/L ~ WW) com- 
ponents, one must typically consider multiple matrix el- 
ement types involving different numbers of W and W' 
fields jfiZ] 

Specifically, wc define 



,(0) Ai.-.A 



ai...a f 



(r)] 



■7p L V [Pl-Pp] 



7 1 --^ tAp 7,» ! 



(5.8) 



where the double angle brackets ((O ■ ■ ■ )) signify the one- 
particle irreducible matrix element of O . amputating the 
external fields. 58 ! 68 The external {W"' } fields are assumed 
to be located far from each other, and from the position 
r of the composite operator. The superscript (0) on the 
left-hand side (LHS) of this equation indicates that this 
is a bare (i.e. not yet renormalized) quantity. 

The basic one-loop process is illustrated in Fig. [3] 
Fig. EJa) depicts the two-field matrix element of the (un- 
symmetrized) operator 



(5.9) 



The vertex QJ4 [Fig.[^b) and Eq. (|5.2|) ] pairwise permutes 
the lower indices of composite operator "legs," as shown 
in Fig. [3jb) . The completely antisymmctrized operator 
defined by Eq. (|2.23p is clearly an eigenoperator at one 
loop, since the sum of all diagrams to this order repre- 
sents a complete symmetrization procedure. This is ex- 
pected to hold to all higher orders in t, because O ^""fl 1 
plays the role of a "highest weight state" in an irreducible 
representation of the full NLerM target manifold symme- 
try group U(2n). 

At one loop, the matrix element defined by Eq. ([57 
is equal to 



-p(O) Ai.-.Ap r^-, oi-a p / \i 
L p 7l ... 7p 1 U p [f} 1 ...fj p \\ r )\ 



1 - 



■h 



21 



Ai ...A p ; Qi ...a p 

V Ti—7pj [Pi--&p] ' 



(5.10) 



where 21 



is the zeroth order amplitude 



Ai ...A p ; ai...a p 
*p 71. ..7 P ; [0!.../3 p ] 

[equal to zero or the pure constant (1/p!) 2 , depending 
upon the matrix clement]^ In order to save writing 
wherever possible, from this place forward we will adopt 
the following shorthand notation: underlined vertices 
(r), operators (O), and tree level matrix elements (21) 
should be understood as possessing the appropriate set 
of indices, and all indices in a given equation are matched 
(in the appropriate order). With these conventions estab- 
lished, Eq. (|5.10|) may be rewritten compactly as 



Lp° } [O p (v) 



pip - 1) 
1 ^ J i 



%■ 



(5.11) 



In Eqs. (IBTOI and ifBTTTp . 

-t f d d k 



h = 



2 
t 

in 



{2n) d |k| 2 + Mo 



1 



hi 



htei 



(5.12) 



Here, we have used Eq. (|5.3p and standard dimen- 
sional rcgularization technology; 7 denotes the Euler- 
Maschcroni constant. 

To renormalizc Eq. (|5.11j) . we insist that^ 



Z- 1 Z#%V\O p {t)]= finite, 



(5.13) 



where Z p is the composite operator rcnormalization, and 

the factor of Z^ 2 compensates for the p (amputated) 
external fields. One then obtains the scaling dimension 

dhi Z D 



dint 



t 

An 



np 



Pip - 1) 



Of 



(5.14) 



We make two observations. First, by evaluating 
Eq. (|5.14[) at the Anderson transition critical point lo- 
cated by t* — 47rv / 2e, with n — > 0, we recover the 
known result s 32 ' 34 for the multifractal spectrum f(p) 
associated to the averaged IPR, as provided above in 
Eqs. (|1.5|) - (|1.6c[) . Second, we have only considered the 
rcnormalization of the fully antisymmctrized operator de- 
fined by Eq. (|2.23p . because this is the most relevant 
in the n — > limit. The fully symmetrized operator 
^P^Pifo '"/3 )' wn i° n * s defined as in Eq. (|2.23p without 
the sgn(P) factor in the summand, also constitutes an 
eigenoperator with scaling dimension 



•Kp (sym) 



t 

-in 



np 



PiP - 1) 



o(t 2 



(5.15) 



Consider the case of n = 1. The fully antisymmctrized 
operator defined by Eq. (|2.23p does not exist for p > 
1, since all W => W fields are scalars in this case. At 
the fixed point located by t* = ine/n with n = 1, the 
symmetrized operator scaling dimension a;* ( - sym - ) = ep(p+ 

l)/2 + O (e 2 ), which is the expected result for spherical 
harmonic composite operators in the U(2)/U(l)xU(l) ~ 
0(3)/0(2) NLaM^l 



C. Two-point function normalization at the MIT 



The set of coefficients {C q q + /} [Eq. OH] defining the 
operator product expansion (OPE) for properly normal- 
ized composite eigenoperators constitute universal num- 
bers characterizing the MIT. The proper (RG scheme- 
dependent) normalization of each eigenoperator is such 
that its two-point autocorrelation function is scheme- 
independent at the critical point in d = 2 + e^ 2 - In this 
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subsection, we derive the normalization of the operators 
defined by Eq. (|2.23|) with respect to their two-point func- 
tions (at large spatial separation), while the OPE is con- 
sidered in the sequel. Since we are interested in critical 
properties, we assume h = in Eq. (|2.7[) throughout the 
following discussion. 

The technical tool for computing operator correlation 
functions at any perturbatively accessible fixed point is 
RG-improved perturbation theory (PT). In the case of 
the non-trivial NLcrM critical point in d = 2 + e, how- 
ever, some technical difficulties arise: For a NLcrM with 
a compact, non-Abelian symmetry, the trivial fixed point 
located at t = is invariably infrared (IR) unstable in 
2Dj£& For any nonzero t, such a 2D model always flows 
under the RG toward a symmetry-restored, thermally- 
disordered "paramagnetic" state. Renormalized pertur- 
bation theory at t <C 1, for composite operator correla- 
tors that are not invariant under the full symmetry group 
of the sigma model target manifold, is typically plagued 
by IR divergences, and hence affected by the specific way 
one regularizes these IR divergences. 

The solutio n 70 i 71 i 72 that we employ in this section is 
to consider only invariant correlation functions. Invari- 
ant correlators are free of IR divergences, and a sensible 
renormalized PT for these objects can be constructed^ 
For the OPE in the next section, we will see that this is 
restriction is unnecessary. 

We stress that, by the same token, all eigenoperators 
at the non-trivial fixed point in d = 2 + e possess well- 
defined critical correlations. Thus the above-described 
calculational impasse, as well as its solution, in fact re- 
flect peculiarities of the e-cxpansion, rather than the 
NLcrM itself (at least for d > 2). 

Within the symmetry-broken phase, V (W^V^) and 

V (\/l n - WW\\Jl n - W^W) fields possess very dif- 
ferent correlation functions: the former constitute Gold- 



stone modes with massless correlations, while the latter 
are gapped longitudinal modes, with massive correlation 
functions for all t < t* . At the critical point t = t* 
for n > 0, symmetry is restored; here, all operators be- 
longing to a given irreducible representation of the target 
manifold symmetry group will possess identical correla- 
tions, provided a group-invariant normalization is chosen 
for these operators. [For the 0(3)/0(2) model, an in- 
variant normalization is that conventionally assigned to 
spherical harmonics, written in terms of tt and a coor- 
dinates.] We will use this fact to determine the two- 
point function normalization of Cp" 1 ^"^ ] [Eq. (|2.23[) ] 
for ?? = {1, 2, . . .}, and then continue the result to n — > 0. 
Consider the following invariant, "non-local" operator, 



fi p (r,r') = ^6 p{m} (r)e; {m} ( 

{m} 



(5.16) 



where O p { m } is a composite operator that is a com- 
ponent of an irreducible representation of the sigma 
model symmetry group. The representation is distin- 
guished by the Casimir parameter p. while the com- 
ponent operators are labeled by a set of "magnetic" 
quantum numbers {m}, e.g. {ot\, . . . , a p , [/?i, . . . , (3 P ]} £ 
{m} for the antisymmetrized LDOS moment operators 
defined by Eq. (|2.23[) . One may use expressions for 
the 6 p { m } in terms of the target manifold coordinates 

(W, W\ Vin - WWi, Vtt - W^W) in order to con- 
struct an explicit expression for the RHS of Eq. (15.161) : 
at the non-trivial critical point in d = 2 + e, however, we 
require only the lowest order expansion for the expecta- 
tion of Eq. (|5.16[) in powers of t* oc e CT (cr = 1 or 1/2 for 
?i G N or 77, — ► 0, respectively). Since an expansion in 
powers of t is equivalent to an expansion in powers of W 
and W\ we make the following ansatz: 



Q„(r,r') ~ 6L fl + ^Tr \w*(t)W(t') - WUv)W(r) + WHt')W{t) - WUr')W(r') 
\ 2n z L 



o(w l 



(5.17) 



The expectation of the assumed form of Eq. (|5 . 1 7[) is free 
from IR divergences. While the numerical coefficient a p 
in this equation is arbitrary, a p \ is not, and can in prin- 
ciple be computed from knowledge of the representation 
theory of the NLcrM symmetry group; instead, we will de- 
termine its value empirically, below. Note that Eq. (|5 . 1 7[) 
is manifestly invariant under U(n)xU(n) subgroup trans- 
formations, W — ► UlWUr, with U\j R U L j R = In- 



We compute the expectation of Eq. (|5 . 1 7[) in position 
spaced using the IR-convergent Green's function at zero 



coupling 



(W^i Pl (r)W a *p 2 (0) - W +ai /j 1 (0)W aa / 5 8 (0))o 



d d k e lk r - 1 



(2ir) d |k| 2 
t 1 



* * (d-2)S d \r\ d ~ 2 



02 01 2lT 



i-Iln(V|r|V) 



(5.18) 



where Sd is the surface area of the sphere in d dimensions, 
and we have used Eqs. (|5.ip and (|5.3|) . 
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Let us define the renormalizcd operator 
[Op] (r,r')^Z- 2 ^( r ' r ') ; 



(5.19) 



where Z p is the renormalization factor obtained via 
Eqs. (|5. 13[) and (|5 . 14[) for the composite operators defined 
by Eq. (|2.23j) ; all operators belonging to a particular irre- 
ducible representation receive the same renormalization 
in a NLctM^I 

Insisting that ([fi p ] (r, r')) is finite (for r ^ r'), we see 
that we must take 



np 



P(P - 1) 



(5.20) 



where y = r r', R = (r + r')/2, a;* is the scaling dimen- 
sion defined by Eqs. ([5^14]) and p^aj) - (fOc|) . and 
is the (universal) OPE coefficient that we seek, expected 
to possess an expansion in powers of t*(e). Eq. (|5.23[) 
will hold as a replacement rule in the limit |y| — ► 0, valid 
inside correlation functions involving arbitrary configu- 
rations of other, spatially remote operators. Note that 
Eq. (|5.23[) relates a product of maximally antisymmet- 
ric operators to a single, maximally antisymmetric op- 
erator; the ellipsis ". . ." on the right-hand side (RHS) 
of this equation represents other, less relevant operators 
that are produced in the 'fusion' process; we will ignore 
the contribution of the latter to functional RG.— 



in Eq. (|5.17[) . We have obtained the lowest order ex- 
pansion coefficient for the group-invariant structure de- 
fined by Eqs. (|5.16[) and (| 5 . 1 7[) without explicitly employ- 
ing group theory, but using only the rcnormalizability of 
the NLerM! [Basic group theoretic knowledge was neces- 
sary to identify the invariant scaling operators defined by 
Eq. (E22J), however.)Zi 

Finally, we set t = t*(e), and then we re-exponentiate 
the expectation of Eq. (|5. 19(1 to obtain, at the non-trivial 
critical point, 



<[fi P ] (r,r')> 



, V 7r^ 2 e 7 ) 



(5.21) 



In order to determine C p + p f 

trary matrix element T p+p , 7 *" \J 

Eq. (f5T23| , as defined by Eq. fet 
and (|5.22p . and employing the compact notation intro- 
duced above and implemented in Eq. (|5.11| . the RHS of 
Eq. (|5.23[) may be written as 



we compute an arbi- 
p+p' r. . . 1 f both sides of 
Using Eqs. ([530)) 



r p °l P ,[O p+pl (R)], 



u p,p' 



I l x i+ x *r— x tj_ i 
y p v p+p 



(5.24) 



where x* is the scaling dimension in Eq. (|5. 14[) . evaluated 
at t = t*. [x* is given explicitly by Eqs. (|1.6a|H1.6cp f° r 
the limit n — > 0, appropriate to the MIT.] 

Eq. (|5.21|) allows us to define the following "renormal- 
ized and normalized" composite operators, which we will 
enclose with the double square brackets [•••]• Referring 
to Eq. (|2.23|) . we designate 



p [01-M 



(r). (5.22) 



where Z\y (Z p+p i) is the field strength (composite oper- 
ator) renormalization factor, and the bare amplitude is 
given by 



r (0) iq 
—p+p' i^-p+p' 



l--(p+p')(p+p'-l)h 
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£±p+p' 



(5.25) 



Eqs. (|5.19[) and (|5.21[) guarantee that the two-point cor- 
relation function between distant operators defined by 
Eq. (|5.22[) is both ultraviolet finite, and independent of 
the renormalization scheme. 



In Eq. (pT25)) . h is the loop integral defined by Eq. i[5"7T2]) . 
where we retain the infrared regularization for now, ho ^ 
0, while 



D. Operator product expansion at the MIT 



X 1 ...X p+p ,;a 1 ...o tp a' L ... a ' p , 
-P+P' ^P+P'fl---7 P+P r, [P 1 .../3 P P' 1 ...P' I ] 



We conclude this section with the construction of the 
operator product expansion (OPE) for the operators in 
Eq. ([5T22]) . At the critical point in d = 2 + e, the OPE is 
expected to take the form 



cp + p' i i 

, ,x;+x;,-x; , I u p+p> lPi...p P P{.. .p',]l W 



(5.23) 



denotes the zeroth order amplitude for the matrix ele- 
ment (a pure number). (See footnote [6^ for details.) 



Similarly, the LHS matrix element of Eq. (|5 . 23[) may 
be written 



p+p' 



xr p °l pl \O p (r)O pl (r% (5.26) 
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r r 



a ( //p ( a 2 //p 2 <V s H3 

T r ^ 




FIG. 4: Basic renormalization process in the OPE. 

where the bare (unrenormalized) amplitude is 
4%,[<2 P (r)O p ,(r')] 

x^t\-^W )) 



p + p 

(5.27a) 



x W ^ ,... W "' P > 0l J r ') 

— 1) other terms 1 
obtained by permutations J 



(5.27b) 



As in Eq. (|5.8p . the external fields {VP^} in these equa- 
tions are assumed to be located far from the vicinity 
of the operator product (r or r') and from each other, 
while the double angle brackets instruct us to take the 
one-particle irreducible amplitude, with external legs 
amputated^ 

Consider the one-loop renormalization of the term 
written explicitly in Eq. (|5.27b[) ; the other (p\p'\ — 1) 
terms implied in this equation will give identical contri- 
butions. The basic renormalization process of an oper- 
ator product is illustrated in Fig. Q] specifically for the 
combination 

W ai ^(r) ® W a ^ 2 (v')W a ^ 3 (r'). 

The vertex QJ 4 [Fig. [2] and Eq. ([5T2]) ] modifies the op- 
erator product in two ways. First, it renormalizes the 
constituent operators, pairwise permuting indices of legs 
both tied to either r or r', as shown in Fig. QJb) [c.f. 



Fig. [3]. Second, QJ4 ties the two operators together 
by pairwise crosspermuting their indices, as depicted in 
Figs. He) and (d). 

Now, algebraically we may express an unsymmetrized 
product of q IP-field matrix elements (W a p) in terms 
of the completely antisymmctrized product, native to 
the most relevant irreducible representation with Casimir 
parameter q of the NLerM symmetry group, plus other 
terms which belong to other (completely symmetric or 
mixed symmetry) representations. In particular, 



W ai Pl ■■■W a "i3 q = w° 



[01 



(5.28) 



where again the square brackets [• ■ ■ ] denote complete 
antisymmetrization. The unity coefficient in front of 
the completely antisymmctrized tensor on the RHS of 
Eq. (|5.28[) follows from the fact that the antisymmetriza- 
tion procedure is a projective operation that kills sym- 
metric or mixed symmetry terms, but leaves the pre- 
antisymmctrized component invariant. 

After the vertex acts upon the unsymmetrized opera- 
tor product displayed explicitly in Eq. (|5.27b[) . giving the 
appropriate factors for the two types of renormalization 
depicted in Fig. 0Jb) an d Figs. BJc),(d), respectively, we 
are free to use Eq. (|5.28|) to replace each resulting un- 
symmetrized, permuted product with the corresponding 
completely antisymmctrized version, up to less relevant 
mixed symmetry or higher gradient terms 2^ The com- 
pletely antisymmetrized product of p + p' factors is just 
the composite operator on the RHS of the OPE, as de- 
fined by Eq. (|5.23p . Therefore, using the Feynman rules 
in Eqs. (|5.ip and (|5.2| . and summing all diagram topolo- 
gies to one loop, Eq. (|5.27a|) may be written as 



p + p' 
P 



1-h 



p{p- I) , p'{p'-l) 



{ -h{y)W] 
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—p+p ' 



(5.29) 



where 2t p+p / denotes the zeroth order matrix element of 



O 



ai ...a p a 1 . 



p+p' [/3i.../3 p /31.../3',p as in Ec L- <1 5 - 25 P> J i is the integral 
defined by Eq. (|FT2]) . and 



-tn 



hnt 



o'o 



(5.30) 



Let us briefly comment upon the origin of the various 
combinatoric factors in Eq. (|5.29|) : The prefactor ( p+ p p ) 
originates from the normalization convention used in 
Eq. (|2~23"1) . The factors of p(p - l)/2 and p'(p' - l)/2 
count the number of inequivalcnt ways leg indices associ- 
ated with either operator O p or Q_ p , (but not both) may 
be permuted, as occured previously in the scaling dimen- 
sion calculation [Eq. (|5.10|) ]. The factor of pp' counts 
the number of inequivalcnt ways one leg index from each 
operator may be intcrpermuted, as in Figs. HJc),(d). 
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Equating the left-hand and right-hand sides of the 
OPE, Eqs. HPHl - flOBD and and expanding ev- 

erything to the lowest non-trivial order in t, we obtain 



cp+p 
°p,p' 



P + p' 



1+pp' [Ji-J 2 (y)]+ln 



Z. 



p+p' 



ZpZp> 



P J — — \Xp+pi ~ Xp — Xpi ) 

xln (Vlyl'e 7 ) 

(5.31) 

We may now take the limit h — ► 0, because the combi- 
nation 2[h - I 2 (y)} [Eqs. (i5TT2"j) and (|Q01) ] gives the IR- 
finite integral evaluated previously in Eq. (|5.18p . above. 
Using Eq. (|5 . 14[) . Eq. (|5.31[) simplifies to the expression 



(jp+p 
p,p' 



P + P 
P 



1 

1 - - 
2 



— x„ 



'p+p' "p 

.2|..|2 -y 



Xp, + t — 
P 47T 



x In (7iy|y|V 

(5.32) 

In general, the RHS of this equation is a UV-fmite, non- 
zero function of the operator separation y. At the non- 
trivial critical point t = t* in d = 2 + e, however, we have 
[from Eq. (|5~14|) ] 



r<p+p' 
°p,p' 



P + p' 
P 



O (t 



* 2^ 



(5.33) 



At the Anderson metal- insulator transition (n — > 0), 
t* = 4ttV27 + O (e). Thus the OPE coefficient for the 
operators with normalization determined by Eq. (|5.22p is 
independent of y/e to order 

The result in Eq. (|5.33|) should be contrasted with a 
similar computation in <f> theory in d = 4 — e: at the 
Wilson-Fisher fixed point, the fusion of two elementary 
renormalized and normalized [0] fields into the mass op- 
erator [c/> 2 ] yields an OPE coefficient that acquires a cor- 
rection at the lowest non-trivial order in the quartic cou- 
pling strength A* oc e + O (e 2 ). (See, e.g., Ref. IzU) As 
in the above NLtrM calculation, the normalization of the 
operators [0] and [0 2 ] is chosen so as to give two-point 
autocorrelation functions independent of the rcnormal- 
ization scheme. 



Finally, tracing over pairs of indices in Eq. (|5 . 23[) allows 
the OPE to be written as 



(jp+p 

[Op] (r) [O pl ] (r') ~ , a . + ff-*V ; ( R ) 

y p p p+p 

where 



a.i — 1 a.p — 1 

= Z- 1 (n^e^ O p (r) 



(5.34) 



(5.35) 



is the renormalized and normalized version of the bare 
operator defined by Eq. (|2.24[) . 

The OPE in Eq. ([534)) . together with the coefficient 

Cjpf given by Eq. (|5.33|) , constitutes the primary tech- 
nical result of this paper. We have succeeded in asso- 
ciating a unique, properly normalized operator [O p ] (r) 
to the p moment of the LDOS, and demonstrated that 
the family of such operators obeys the OPE set forth 
in Eq. (|1.11[) in the Introduction. This is the necessary 
input to the functional RG scheme used in Sec. Illll to ex- 
tract the typical multifractal spectrum r(q) in Eq. (|3.24l) . 
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addition formula 

= 2_. {hh; lm\lih; mim2) c\ lth Yi tm (n + ,n-,a), 

l ,m 

where (£1/2; lm\l\_l2\ m\m?) is a Clebsch-Gordan coefficient, 
and cj j is a normalization constant that is independent 
of {m, mi, 7712}. As an example, one has 



(7T + +7T_) = — 
A2 



*2,2 - \l -Y 2 ,0+Y 2 ,-i 



1 4 V 

+ T— o ^0,0, 
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where the coefficients {A2,Ao} are determined by 
Eq. dUTHJl. 

Note, however, that for the 0(3)/0(2) model, higher I 
operators possess less relevant scaling dimensions at the 
non-trivial critical point in d — 2 + e.— The existence of 
operators with arbitrarily negative scaling dimensions in 
the unitary class NLerM emerges due to the replica n — » 
limit. 

Because of the symmetry restoration linking the l n' 

(W,W*) and V (VL - WW*, VL - WW) NLaM 
fields at the non-trivial critical point in d = 2 + e (for 
fixed n£ {1,2, ...}), many other eigenoperators exist that 
share the same "maximally" relevant eigenvalues as those 
defined by Eqs. (|2.23|l and (|2.24[1 , Our choice to focus upon 
the latter operator family is one of convenience. 
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It seems that a Wilsonian RG approach, such as that im- 
plemented in Ref. l44l . does not produce as easily the fea- 
tures that we seek (such as, e.g. the lowest order fusion 
process 010 01 — » O2, as represented by a term cx Y\ 
in the one-loop beta function for Y%). Rather than using 
the FRG within the Wilsonian scheme, we deploy here the 
field theory approach explicated in Sec. [V] 
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Strictly speaking, the (bare, unrenormalized) matrix ele- 
ment in Eq. (|5.8p is a mixed position and momentum space 
object, the precise definition of which is given by the cor- 
relation function 
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where G(k) is the "one-particle" propagator, which has 
the amplitude given by Eq. (|5.1|l at tree level. Only com- 
pletely connected diagrams are counted in the evaluation 
of the right-hand side of Eq. (|5.36|l . The external momenta 
{ki, . . . , k p } can be taken to vanish in this equation. 
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ten as a sum of 2p-fold products of Kronecker delta 
functions. In the case of the operator 2 °j^^](r) = 
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' 3 As discussed in Ref. |72i . the basic problem is that non- 
invariant correlation functions suffer ultraviolet (UV) di- 
vergences for d > 2, and infrared (IR) divergences for 
d < 2 — there is no finite window of dimensionality free 
from both UV and IR problems. Within dimensional regu- 
larization, non-covariant correlators cannot be regularized 
via analytical continuation in d, so renormalized PT does 
not exist for these objects. By contrast, invariant correla- 
tors can be shown to be IR and UV finite for < d < 2; 
in this case, dimensional regularization works in the usual 
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the consequence of the Ward identity. [See also the discus- sity Press, New York, 1992). 

sion following Eqs. (fT3]) - (fOb)) .] 
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